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The detection of the theoretically expected dark matter is central to particle physics 
and cosmology. Current fashionable supersymmetric models provide a natural dark 
matter candidate which is the lightest supersymmetric particle (LSP). Such mod- 
els combined with fairly well understood physics like the quark substructure of 
the nucleon and the nuclear structure (form factor and/or spin response function), 
permit the evaluation of the event rate for LSP-nucleus elastic scattering. The 
thus obtained event rates are, however, very low or even undetectable. So it is 
imperative to exploit the modulation effect, i.e. the dependence of the event rate 
on the earth's annual motion. Also it is useful to consider the directional rate, i.e 
its dependence on the direction of the recoiling nucleus. In this paper we study 
such a modulation effect both in non directional and directional experiments. We 
calculate both the differential and the total rates using both isothermal, symmetric 
as well as only axially asymmetric, and non isothermal, due to caustic rings, veloc- 
ity distributions. We find that in the symmetric case the modulation amplitude is 
small. The same is true for the case of caustic rings. The inclusion of asymmetry, 
with a realistic enhanced velocity dispersion in the galactocentric direction, yields 
an enhanced modulation effect, especially in directional experiments. 

1 Introduction 

In recent years the consideration of exotic dark matter has become necessary in 
order to close the Universe a. Furthermore in in order to understand the large 
scale structure of the universe it has become necessary to consider matter 
made up of particles which were non-relativistic at the time of freeze out. 
This is the cold dark matter^omponent (CDM). The COBE data □ suggest 
that CDM is at least 60% el. On the other hand during the last few years 
evidence has appeared from two different teams, IheJJigh-z Supernova Search 
Team □ and the Supernova Cosmology Project tl ' □ which suggests that the 
Universe may be dominated by the cosmological constant A. As a matter 
of fact recent data the situation can be adequately described by a baryonic 
component Qb = 0.1 along with the exotic components Qcdm = 0.3 and 
51a = 0.6. In another analysis Turner 13 gives Q m = flcDM + Ufl n = 0.4. 
Since the non exotic component cannot exceed 40% of the CDM El, there 
is room for the exotic WIMP's (Weakly Interacting Massive Particles). In 
fact the DAMA experiment □ has claimed the observation of one signal in 
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direct detection of a WIMP, which with—better statistics has subsequently 
been interpreted as a modulation signal tl 

The above developments are in line with particle physics considerations. 
Thus, in the currently favoured supersymmetric (SUSY) extensions of the 
standard model, the most natural WIMP candidate is the LSP, i.e. the lightest 
supersymmetric particle. In the most favoured scenarios the LSP can be 
simply described as a Majorana fermion, aJinear_cpmbination of the neutral 
components of the gauginos and Higgsinos ldOE3o. 

Since this particle is expected to be very massive, m x > 30GeU, and 
extremely non relativiatic. with average kinetic energy T < lOOKeV, it can 
be directly detected ESIlj mainly via the recoiling of a nucleus (A,Z) in the 
elastic scattering process: 

X + (A,Z) - X + (A,Z)* (1) 

(x denotes the LSP). In order to compute the event rate one needs the fol- 
lowing ingredients: 

1) An effective Lagrangian at the elementary particle (quark) leveLob- 
tained in the framework of supersymmetry as described , e.g., in Refs. uej. 

2) A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of the 
nucleon in quarks other than u and d. This is pajdiciilarly true for the scalar 
couplings as well as the isoscalar axial coupling E£rEL 

3) Compute the relevant nuclear matrix elements ErE3 using as reliable as 
possible many body nuclear wave functions. By putting as accurate nuclear 
physics input as possible, one will be able to constrain the SUSY parameters 
as much as possible. The situation is a bit simpler in the case of the scalar 
coupling, in which case one only needs the nuclear form factor. 

Since the obtained rates are very low, one would like to be able to exploit 
the modulation of the event rates due to the earth's revolution around the 
sun E§Ea~E3 Jb^this end one adopts a folding procedure assuming some 
distribution EJj23c3 of velocities for the LSP. One also would like to know 
the directional rates, by observing the nucleus in a certain direction, which 
correlate with the motion of the sun around the center of the galaxy and the 
motion of the Earth BE3. 

The calculation of this cross section has become pretty standard. One 
starts with representative ppiit in the restricted SUSY parameter space as 
described in the literature E3a. We will adopt a phenomenogical procedure 
taking universal soft SUSY breaking terms at Mqut, i-e., a common mass for 
all scalar fields mo, a common gaugino mass Afjya and a common trilinear 
scalar coupling Aq, which we put equal to zero (we will discuss later the 
influence of non-zero Ao's). Our effective theory below Mqut then depends 
on the parameters E3: 

mo, M1/2, Mo, ac Mgut, h t , , h b , , h T , tan/3 , 
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where aa — g^/iTr (gG being the GUT gauge coupling constant) and h t ,hb, h T 
are respectively the top, bottom and tau Yukawa coupling const ant Sj-ajt Mqut- 
The values of ac and Mqut are obtained as described in Ref.LS. For a 
specified value of tan/3 at Ms, we determine h t at Mqut by fixing the top 
quark mass at the center of its experimental range, mt{mt) — 166GeV. The 
value of h T at Mqut is fixed by using the running tau lepton mass at mz, 
m T (mz) — 1.746GeV. The value of hb at Mqjjt used is such that: 

= 2.90 ± 0.14 GeV. 

after including the SUSY threshold correction. The SUSY parameter space is 
subject to the following constraints: 

1. ) The LSP relic abundance will satisfy the cosmological constrain: 

0.09 < n LSP h 2 < 0.22 (2) 

2. ) The Higgs bound obtained from recent CDF @ and LEP2 @, 
i.e. m h > 113 GeV. 

3. ) We will limit ourselves to LSP-nucleon cross sections for the scalar cou- 
pling, which gives detectable rates 

4 x 10- 7 pb < a™t™ < 2 x lO" 5 pb (3) 

We should remember that the event rate does not depend only on the nucleon 
cross section, but on other parameters also, mainly on the LSP mass and the 
nucleus used in target. The condition on the nucleon cross section imposes 
severe constraints on the acceptable parameter space. In particular in our 
model it restricts tan(3 to values tan[3 ~ jaQ. We will not elaborate further on 
this point, since it has already appeared 113. 



2 Expressions for Extracting of the Nucleon Cross Section 
from the Data. 



The effective_Lagrangian describing the LSP-nucleus cross section can be cast 
in the form E3 

G 

L eff = — ^|{(xi7 A 75Xi)^A + (xiXi)J} (4) 

where 

Jx = Nixifv + fvr 3 + f° Al 5 + ffo 5 r 3 )N , J — N(fs + f}r 3 )N (5) 

We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, xi7 A Xi — (identically). 
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With, the,, above ingredients the differential cross section can be cast in 
the form BEEil 

dii v 
da{u,v) = — [(Eg + t v -2) F 2 (u) + t spin F 11 (u)] (6) 

Ss - ^Miv) } {A [ifs ~ fs ^A~ ] ] ~ a ^ oA ( Miv) } (7) 

yi spin j- /q\ 

^spin — &p ^,0 ^spin 

3(1 + ^a)2 fX F n (u) fX Fn(u) 

J A 

(9) 

t v = a v v x0 C v (10) 

_ ( t i r (A)/p r (N)f A 2 h fv A-2Z 2 v 2 1 2^+1 (2u) 

W (1 + 4 )2 A [i /° A c (2rf(l + 7 7) 2 (« 2 > J 

(11) 

cr* x o = proton cross-section, i = S 1 , spin, 1/ given by: 

<7p x0 = Co (is) 2 ) 2 (scalar) , (the isovector scalar is negligible, i.e. 

" spm ' 3 (/° + /A) 2 (^f) 2 (spin) , < x0 = a (f v + ftf (^f>) 2 



p,x" 
(vector) 

where m^r is the nucleon mass, rj = m x /mNA, and n r {A) is the LSP-nucleus 
reduced mass, /z r (iV) is the LSP-nucleon reduced mass and 

°"o = ^-(G F m N ) 2 ~ 0.77 x 10~ 38 cm 2 (12) 

27T 

U = g 2 6 2 /2 , Q = Q « , ^ = ^7^2 (I 3 ) 

where b is (the harmonic oscillator) size parameter, q is the momentum trans- 
fer to the nucleus, and Q is the energy transfer to the nucleus. 
In the above expressions F(u) is the nuclear form factor and 

are the spin form factors (p, p are isospin indices) Both form factors are 
normalized to one at u — 0. f^o (^1) are the static isoscalar (isovector) spin 
matrix elements. 
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The non-directional event rate is given by: 

n ^ dN p(0) m 

R = Rnon-dir = —rr = a <T{U,V)\V\ (15) 

at m x ArriN 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the 
detector mass The differential non-directional rate can be written as 

dR = dR non - dir = m da(u,v)\v\ (16) 

where da(u,v) was given above. 1—11— i 

The directional differential rate lilljHa in the direction e is given by : 

dR d ir = ^-^—v.eH{v.e) A. da(u,v) (17) 
m x ArriN 27r 

where H the Heaviside step function. The factor of 1/27T is introduced, since 
the differential cross section of the last equation is the same with that entering 
the non-directional rate, i.e. after an integration over the azimuthal angle 
around the nuclear momentum has been performed. In other words, crudely 
speaking, l/(27r) is the suppression factor we expect in the directional rate 
compared to the usual one. The precise suppression factor depends, of course, 
on the direction of observation. In spite of their very interesting experimental 
signatures, we will not be concerned here with directional rates. The mean 
value of the non-directional event rate of Eq. (jlj) , is obtained by convoluting 
the above expressions with the LSP velocity distribution f(v, ve) with respect 
to the Earth, i.e. is given by: 

/dR\ p(0) m f ^ i ^!h^l d 3 v ( lg ) 
\ du I m x ArriN J ' du 

The above expression can be more conveniently written as 

ldR\ p(0) m rj-gr dE, 



where 

After performing the needed integrations over the velocity distribution, 
to first order in the Earth's velocity, and over the energy transfer u the last 
expression takes the form 

R = R t [1 + h(a, Q min )cosa)] (21) 

where a is the phase of the Earth (a = around June 2nd) and Q m in is 
the energy transfer cutoff imposed by the detector. In the above expressions 
R is the rate obtained in the conventional approach L_3 by neglecting the 
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folding with the LSP velocity and the momentum transfer dependence of the 
differential cross section, i.e. by 




tains all the parameters of the SUSY models. The modulation is described 
by the parameter h . Once the rate is known and the parameters t and h, 
which depend only on the LSP mass, the nuclear form factor and the veloc- 
ity distribution the nucleon cross section can be extracted and compared to 
experiment. 

The total directional event rates can be obtained in a similar fashion by 
by integrating Eq. (|l^) with respect to the velocity as well as the energy 
transfer u. We find 



We remind that the z-axis is in the direction of the sun's motion, the y-axis is 
perpendicular to the plane of the galaxy and the x-axis is in the galactocentric 
direction. The effect of folding with LSP velocity on the total rate is taken 
into account via the quantity t°, which depends on the LSP mass. All other 
SUSY parameters have been absorbed in R. We see that the modulation of 
the directional total event rate can be described in terms of three parameters 
hi, 1=1,2,3. In the special case of A = we essentially have one parameter, 
namely hi, since then we have = 0.117 and /13 = 0.135. Given the functions 
hi(a, Qmin) one can plot the the expression in Eq. ( p3| ) as a function of the 
phase of the earth a. 

3 The Scalar Contribution- The Role of the Heavy Quarks 

The coherent scattering can be mediated via the the neutral intermediate 
Higgs particles (h and H), which survive as physical particles. It can also be 
mediated via s-quarks, via the mixing of the isodoublet and isosinlet s-quarks 
of the same charge. In our model we find that the Higgs contribution becomes 
dominant and, as a matter of fact the heavy Higgs H is more important 
(the Higgs particle A couples in a pseudoscalar way, which does_not lead to 
coherence). It is well known that all quark flavors contribute Eu, since the 
relevant couplings are proportional to the quark masses. One encounters in 
the nucleon not only the usual sea quarks (uu, dd and ss) but the heavier 
quarks c, 6, £_which couple to the nucleon via two gluon exchange, see e.g. 



Rd, 



R[(t°/4w) |(1 + h\{a, Q min )cosa)e z .e 

h%{a, Qmin) cosae v .e + h 3 (a, Q m in) sinae x .e\ 



(23) 
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As a result one obtains an effective scalar Higgs-nucleon coupling by using 
effective quark masses as follows 

m u — > f u rriN, m d ^f d m N . m s -> f s m N 
rtiQ — > fQiripf, (heavy quarks c, b, i) 
where m^r is the nucleon mass. The isovector contribution is now negligible. 
The parameters f q , q — u,d,s can be obtained by chiral symmetry breaking 
terms in-celation to phase shift and dispersion analysis. Following Cheng and 
Cheng 123 we obtain: 

/„ = 0.021, f d = 0.037, f a = 0.140 (model B) 
f u = 0.023, f d = 0.034, f s = 0.400 (model C) 
We see that in both models the s-quark is dominant. Then to leading order 
via quark loops and gluon exchange with the nucleon one finds: 

f Q = 2/27(1- £,/«) 

This yields: 

f Q = 0.060 (model B), f Q = 0.040 (model C) 
There is _a_ correction to the above parameters coming from looys. involving 
s-quarks 113 and due to QCD effects. Thus for large tan/3 we find lid: 
f c = 0.060 x 1.068 = 0.064, f t = 0.060 x 2.048 = 0.123, f b = 

0.060 x 1.174 = 0.070 (model B) 
f c = 0.040 x 1.068 = 0.043, f t = 0.040 x 2.048 = 0.082, f b = 
0.040 x 1.174 = 0.047 (model B) 
For a more detailed discussion we refer the reader to Refs 

4 Results and Discussion 

The three basic ingredients of our calculation were the input SUSY parameters 
(see sect. 1), a quark model for the nucleon (see sect. 3) and the velocity 
distribution combined with the structure of the nuclei involved (see sect. 2). 
we will focus our attention on the coherent scattering and present results for 
the popular target 127 1. We have utilized two nucleon models indicated by B 
and C which take into account the presence of heavy quarks in the nucleon. 
We also considered energy cut offs imposed by the detector, by considering 
two typical cases Q m in = 10, 20 KeV. The thus obtained results for the 
unmodulated total non directional event rates Rt in the case of the symmetric 
isothermal model for a typical SUSY parameter choice 113 are shown in Fig. 
|]. Special attention was paid to the the directional rate and its modulation 
due to the annual motion of the earth in the case of isothermal models-^ The 
case of non isothermal models, e.g. caustic rings, is more complicated c3 and 
it will not be further discussed here. As expected, the parameter to, which 
contains the effect of the nuclear form factor and the LSP velocity dependence, 
decreases as the reduced mass increases. 

We will concentrate in the case of isothermal models we will limit and 
restrict ourselves to the discussion of the directional rates. In the special case 
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Figure 1. The Total detection rate per (kg — target)yr vs the LSP mass in GeV for a typical 
solution in our parameter space in the case of 127 / corresponding to model B (thick line) 
and Model C (fine line). For the definitions see text. 



of the direction of observation being close to the coordinate axes the rate is 
described in terms of t and the three quantities hi, i = 1, 2, 3 (see Eq. (gl)). 
As expected, the parameter to, decreases as the reduced mass increases. 
These are shown in tables 1-2 for various values of Q m in and A. ft. 2 and /13 are 
constant, 0.117 and 0.135 respectively, in the symmetric case. On the other 
hand hi, hi and /13 substantially increase in the presence of asyjipactry. For 
the differential rate the reader is referred to our previous work qEj. 

5 Conclusions 



In the present paper we have discussed the parameters, which describe, the 
event rates for direct detection of SUSY dark matter. It is well known EH the 
event rates are quite small and only in a small segment of the allowed pa- 
rameter space they are above the present experimental goals. One, therefore, 
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Table 1. The quantities t° and hi for A = in the case of the target 53/ 127 for various 
LSP masses and Q m in in KcV (for definitions see text). Only the scalar contribution 
is considered. Note that in this case h 2 and I13 are constants equal to 0.117 and 0.135 
respectively. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


1.960 


1.355 


0.886 


0.552 


0.442 


0.360 


0.212 


hi 


0.0 


0.059 


0.048 


0.037 


0.029 


0.027 


0.025 


0.023 


t° 


10. 


0.000 


0.365 


0.383 


0.280 


0.233 


0.194 


0.119 


hi 


10. 


0.000 


0.086 


0.054 


0.038 


0.033 


0.030 


0.025 


t° 


20. 


0.000 


0.080 


0.153 


0.136 


0.11 


0.102 


0.065 


hi 


20. 


0.000 


0.123 


0.073 


0.048 


0.041 


0.036 


0.028 



Table 2. The same as in the previous, but for the value of the asymmetry parameter A = 1.0. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


2.429 


1.825 


1.290 


0.837 


0.678 


0.554 


0.330 


hi 


0.0 


0.192 


0.182 


0.170 


0.159 


0.156 


0.154 


0.150 


h 2 


0.0 


0.146 


0.144 


0.141 


0.139 


0.139 


0.138 


0.138 


h :i 


0.0 


0.232 


0.222 


0.211 


0.204 


0.202 


0.200 


0.198 


t° 


10. 


0.000 


0.354 


0.502 


0.410 


0.349 


0.295 


0.184 


hi 


10. 


0.000 


0.241 


0.197 


0.174 


0.167 


0.162 


0.154 


h 2 


10. 


0.000 


0.157 


0.146 


0.142 


0.140 


0.140 


0.138 


h :i 


10. 


0.000 


0.273 


0.231 


0.213 


0.208 


0.205 


0.200 


t° 


20. 


0.000 


0.047 


0.169 


0.186 


0.170 


0.150 


0.100 


hi 


20. 


0.000 


0.297 


0.226 


0.190 


0.179 


0.172 


0.159 


h 2 


20. 


0.000 


0.177 


0.153 


0.144 


0.142 


0.141 


0.139 


h :i 


20. 


0.000 


0.349 


0.256 


0.224 


0.216 


0.211 


0.203 



is looking for characteristic signatures, which will aid the experimentalists 
in reducing background. These arc two: a) The non directional event rates, 
which arc correlated with the motion of the Earth (modulation effect) and 
b) the directional event rates, which are correlated with both the velocity of 
the sun, around the center of the Galaxy, and the velocity of the Earth. We 
separated our discussion into two parts. The first deals with the elementary 
aspects, SUSY parameters and nucleon structure, and is given in terms of the 
nucleon cross-section. The second deals with nucleon and nuclear effects. In 
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the second step we also studied the dependence of the rates on the energy cut 
off imposed by the detector. 

Even though the actual velocity dependence may._arise out of a combina- 
tion of isothermal and non isothermal contributions E3, in the present paper 
we focused on isothermal models. Detectable rates are possible with some 
choices of the input parameters. A typical graph for the total unmodulated 
rate is shown Fig. [I]. We will concentrate more on the directional rates, which 
arc described in terms of the parameters to, hi, hi and h%. To simplify matters 
these parameters are given in Tables 1-2 for directions of observation close to 
the three axes x, y, z. We see that the unmodulated rate scales by the cos6 s , 
with S being the angle between the direction of observation and the veloc- 
ity of the sun, The reduction factor of the total directional rate, along the 
sun's direction of motion, compared to the total non directional rate depends, 
of course, on -the nuclear parameters, the reduced mass and the asymmetry 
parameter A c3. It is given by the parameter f re d = to/ (An to) = k/(2tt). 
We find that k is around 0.6 for no asymmetry and around 0.7 for maximum 
asymmetry (A = 1.0). In other words it is not very different from the naively 
expected f re d = l/(27r), i.e. k = 1. The modulation of the directional rate 
increases with the asymmetry parameter A and it depends, of course, on the 
direction of observation. For Q m ;„ = it can reach values up to 23%. Values 
up to 35% are possible for large values of Q m in, but they occur at the expense 
of the total number of counts. 

This work was supported by the European Union under the contracts 
RTN No HPRN-CT-2000-00148 and TMR No. ERBFMRX-CT96-0090 and 
UENEA 95 of the Greek Secretariat for Research. 
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